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ABSTRACT 

Ostwald-ripening theory is deduced and discussed starting from the 
fundamental principles such as Xsing model concept, Mayer cluster expansion, 
Langer condensation point theory, Ginzburg-Landau free energy, Stillinger 
cutoff-pair potential, LSW-theorv and MLSV-theory. Mathematical intricacies 
are reduced to an understandable version. Comparison of selected works, from 
1949 to 1 98- . on solution of diffusion equation with and without sink/source 
term(s) is presented. Ka'niveit's 1980 work and Marqusee-Ross ' 198- worn are 
more emphasized. Odijk and Lekkerkerker ' s 1985 work on rodlike macromolecules 
is introduced in order to stimulate interested investigators . 


ORIGINAL PAGE IS 
OF POOR QUALITY 


CONTENTS 

I. INTRODUCTION 1 

1.1. Theory of Ostwald-ripening , pioneers' work and development 1 

1.2. Phase separation, its hydrodynamic/chemical-kinetic analogy 3 

1.3. Theories of crystal growth, LSW-theory and MSLW-theory 5 

I. A. Time rate of changing particle radius - growth law 7 

1.5. Works by Marqusee-Ross and by Kahlweit 7 

1.6. The main goal of this work 8 

(Figures 1 & 2), see Figure Captions. 10 

II. FUNDAMENTAL CONCEPTS 11 

11. 1. Diffusion equations 11 

11.2. LSW-theory and Kahlweit's deduction 13 

11.3. Sink density and competition of sinks in precipitation 14 

11. 4. Particle size plays an important role in precipitate competition 15 

11.5. Mass transfer in crystal growth 17 

11.6. Entropy production and Onsager reciprocal relations 19 

11.7. Thermodynamic argument or, cry'stai growth 2i 

111. ASYMPTOTIC DISTRIBUTION FOR PARTICLE SIZES 24 

III. 1 . Fundamental equation and reduced variables 24 

111. 2. Decay function in supersaturation 25 

111. 3. Power expansion and scaled distribution function 26 

111. 4. Scaled growth rate as function of scaled radius 29 

111. 5. Scaled distribution in zero volume fraction of sinks 30 


ORIGINAL PAGE fS 
OF POOR QUALITY 


III. 6. Scaled distribution in nonzero volume fraction of sinks 32 


Table 1. uKz, (X^)^z C^-z 

x=(2+X) * , vs z in different A- and 34 

IV. DISCUSSION AND CONCLUDING REHARKS 35 

IV. 1. Concentration gradient as driving force in diffusion 35 

IV.2. Solution of Ficks first law, comparison between selected works 35 

IV. 3. Approximate solution for complicated process, example 36 

IV. 4. Experimental efforts following the theories 39 

IV. 5. Application of theory to alloys 40 

IV. 6. System containing nonspherical sinks 42 

Table 2. Comparison of selected works on solution of 

d C/2 t=Dv C+A , with A=0 and or A^O 44 

(Figures 3 to 9), see Figure Captions. 44" 

APPENDIX A. Evaluation of the growth law 45 

APPENDIX B. Continuity equation and its Laplace transform 46 

APPENDIX C. Muitiscattering representation and its Fourier transform 53 

APPENDIX D. Solution of F let's 2nd lav, example 59 

REFERENCES 62 

. FIGURE CAPTIONS 69 

(9 figures followed.; 


it 



ORIGINAL PAGE IS 
OF POOR QUALITY 


1. INTRODUCTION 


1.1. Theory of Ostwald-ripening , pioneers" work and development 

Theories for liquid or solid binary mixture-nucleation, coagulation and 
droplet growth have been discussed extensively by Binder and Stauffer [1]. 
The'works of pioneers such as Ising model concept [2], Mayer cluster 
expansion [3], and the Langer-condensation point theory 14], lead to the 
analytical continuation of the free energy, beyond the experimentally 
unobservable singularity [3], by conventional droplet model. 

Ising-spin system in a magnetic field and a lattice gas are proved to 
be mathematically equivalent [2]. Mayer cluster expansion considers that 
the negative sum of the magnetic field and free energy, -(H + F), which is 
equivalent to pressure of the gas, can be expanded into cluster integrals 

N 

_ 1 

b. as p=kT b.z . In which k is Boltzman constant, z is the 
active number density, T is the absolute temperature, and b^ the cluster 
integrals defined over the coordinates of i molecules in volume V. The 
cluster integral depends on the potential energy and coordinates of the i 
molecules [2]. However, the statistical theory of condensation phenomena 
aevelcpec by Mayer and his coworkers and subsequently by Kanr. ant trier-been . 
and born and Fuchs [5j, with its mathematical intricacy, was hardly applies 
to experimental works. 

Accenting to hanger's condensation theory [4;, the stage of degree of 

— i 

evolution is expressed most naturally in terms of the wave number k. At 

k the structure function S(k,t) is a maximum. The structure function 
m 

S(k,t) is defined as the Fourier transform of the spatial correlation function 
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G(r,t). In principle the structure function can be obtained experimentally 
by light scattering. 
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Immediately following the pioneers' works, Frenkel [6] proposed in 
1939 a simplified format of mathematical derivation, which emphasizes that 
the importance of a droplet radius, therefore the velocity of condensation 
of a oversaturated vapor and the growth of the droplet can be deduced by the 
number of molecules (or atoms for monoatomic species) in cluster, since the 
number of molecules in cluster is related to the average radius of the 
clusters each containing the same number of molecules or atoms. 

Hill [7a, b] applied the cluster-theory in 1935 to imperfect gases and 
protein solutions. Greenwood [7c] analyzed theoretically the diffusion flow, 
the smaller to the larger particles for the growth of dispersed precipitates 
in solutions. Perhaps all these works were related to Bigelow and Trimble's 
extensive experimental work [ 7 d ] on vapor pressure and particle size. They 
repeated Ostwaid's experiments [7e] on sulfur, since Ostvald found that, the 
vapor pressure difference between two droplets of unequal size. Sulfur was 
the species Ostvald suggested other ones to use for verifying his finding 
[7d,e]. Finally Greenwood .'7cj tried to verify the particle size and vapor 
pressure relationship, anc tns growth rate, by use of the well known 
Thomson-Freundlicr. equation .'7f ; . All these lead to the name Ostweld-ageir.g 
[9bj and or Ostwalc-ripeninr [lfe,2-a,bj. 

Greenwood [7c; derived the formulas for the rate of change of particle 
size and performed calcuiati onss cr growth rate for highly dispersed particles. 
He applied the theory to experimental results of growth in uranium-lead and 
uranium-sodium slurries. 
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Stillinger (8] proposed in 1963 a concept of "physical cluster integral" 
by overlap of spheres drawn each molecule, which are not regarded as statisti- 
cally independent, but their mutual exclusion property gives rise to free 
energy of formation, work for cavities necessary to contain those clusters. 
Stillinger's work gives a connection between these free energies and a 
suitable generalization of the Mayer's irreducible cluster sum. Above all, 
the most important points are: 

1. Hill included the external forces (such as gravitational field) 
upon the system considered. 

2. Stillinger [8] based on Frenkel-Bank theory of association 
equilibrium intensively and introduced the idea cutoff-potential, 
V(r), which vanishes identically beyund r = b, where b/2 is the 

molecular radius (see Fig. 1). 

The kinetics of precipitation from supersaturated solutions, known as 
Ostwald ripening [ 7 d , e ] , were examined in detail by Lifshitz and Slvozov 
[9a], and by Wagner [9b] independently and conventionally termed LSW-theorv. 
Formulas for the asymptotic grain size distribution, for the number of grains 
per unit volume, and for the supersaturation as a function of time, were all 
established by the LSW-theorv. 

1.2. Phase separ ation . its- h y c r o d y n am i-c / c h e e. 1 n a 1 - k i n e t i c ana _ogy 

A thermodynamically unstable system separates into its components is 
termed spinodal decomposition (10a,bj. Phase separation is clearly c. 
nonequilibrium and nonlinear phynomenon. For this reason, Siggia [10a. used 
models of fluid systems that typically would manifest themselves ir. me so- 
called later stage of spinodal decomposition. Prior to Siggia's work, Langer 
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and hi6 coworkers [ 10b ] performed series of calculations with the theory of 
spinodal decomposition through reagion in which it i6 difficult to distinguish 
from nucleation and growth. The experimental part concerning this topic was 
carried out by Cahn and Hilliard in metallurgical systems [10c,d]. Light 
scattering experiments were carried out by Huang, Goldburg, Bjerkaas, Wong, 
Knobler and Chou I lOe-h] , in order to verfy Siggia's theoretical argument 
[10a] and Langer, Bar-on, and Miller (LBM) theory [10b], More detail can 
be found in Section IV. 


The interaction between molecules or between ions, which leads to 

cluster formation and eventually to the evolution of crystals, has been 

likened to a chemical reaction so that the chemical kinetic theories can 

apply [11a]. Similar to chemical reaction, the activation energy the 

barrier to nucleation must be overcome before crystallization can occur, 

60 that the intermeadiate stage and the last stage of precipitation can 

follow. One result of the energy barrier is the necessary creation of some 

degree of supersaturation before spontaneous crystallization will occur 

[llaii]. The cluster formation can be dealt with chemical equilibrium 

expressions [22c], for instance, £=A. /A C for nA=A and $ = A /A A. 

n n n — k n-t-k n k 

for A +A. =A ... However, in addins: or subtracting monome(s) to or from a 
n k n-i-k ~ ~ 

growing aggregate , thousands of steps forward anc backward may be needed 
to consider . 


Siggia [ 1 0 a j discussed in detail tne influence of hydrodynamic interac- 
tion on the coarsening rate of a mast of cropiets combining througr. diffusive 
coalescence. He emphasized that tne "volume fraction" of precipitate actual! v 
produced in most off-critical quencn experiments probably favors direct 
coalescence (the zero volume fraction was taken as in LSW - theory). 
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According to Landau's theory the free energy can be expanded in a power 
series of the order parameter ^ and it6 derivative. By minimizing the 
free energy one obtains the most probable value of ^ . Landau made 
assumption that the most probable value of ^ coincides with the mean value 
^ . Ginzburg modified Landau's theory to a more realistical form which 
bears the name Ginzburg-Landau free energy [10a]. 

1.3. . Theories of crystal growth, LSW-theory and MLSW-theory 

The theories of crystal growth have been divided into two categories: 

1. Those describing the rate-limiting process in terms of transport to 
the crystal surface are called diffusion controlled (or limited) 
growth (dl-cases). 

2. Those relating the rate control to processes occuring in the 
interface region (i.e., precipitation of crystals from a stirred 
liquid solution) are called interface controlled (or kinetic) 
growth (ik-cases). 

Diffusion controlled growth has been a subject of continued imerest for 
decades [llj, starting probably with the old Noyes-Vhitney [lib], Nernst 
[lie.; equations and developing to the sophisticated approaches by Frisch 
[ 1 2 j , Han. [13], finally by Lifshitz and Slyozov and Vagner [9]. 

The theory of dif f usion-contrcilec particle coarsening developed by 
LSV was modified ny Ardeli jin], to take into account the volume fraction, 

, of the Dtecinitate [15i. The characteristic distance chosen for 
computational purpose was related tc the mean free path between a particle 
and its nearest neighbor. Also spherical symmetry was assumed. It was named 
the modified LSV theory (MLSV-theory) by Ardell [14], The MLSW-theory 


5 



ORIGINAL PAGE* AS 
OF POOR QUALUn? 

1/3 

predicts that while the basic t kinetic of the LSW theory is maintained, 
the coarsening rate increases with increasing volume fraction even at very 
small volumes of [14-16, 21b). 

At zero volume fraction ( =0) the MLSW-theory is equivalent to the 
LSW-theory. The MLSW theory was justified for several alloys such as with 
the available data in the coarsening rates of A-precitates in A-B alloys 
such as Co in Cu-Co alloy. Ardell discussed [14) the works by Sarian and 
Weart and Asimow [16] as they proposed that the mean particle diameter d 

3 V 

and number of particles per unit volume N y are functioin of temperature T, 
volume fraction cj> , and time t, 

d av = f(T, <#> ) t 1/3 , N y = f(T, 4> ) t" 1 

Ardell also emphasized that the Dirichlet region for a 2-dimensional 
distribution of circular particles (see Fig. 2), the concentration of solute 
in the matrix must approach 

C(r) r=r . = C(r') = C' 

at a distance from the ith particle, where 

r' = r, -*■ L/2, L = mean free path 

Feiderhof and Deutch [ 1 7 ] dealt with the concentration aepender.ee of tne 
rate of dif f usior.-contrcliec reactions, particularly in steady state situation 
for a random distribution of the sinks (precipitate particles;. The sir.x 
concentcrstior. is a function of rate constant k, the size ccstrioutior. 
function F . Tne rate constant is predicted to increase with concentration of 
sinks and the dependence on concentrac ionis shown to be nonanalytic. 
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1.4. Tine rate of changing particle radius-growth law 


Brailsford and Wynblatt [18] derived the spatially homogeneous rate 

theory model, to describe the time rate of change of radius of a spherical 

particle embedded in a conf gurational random array of particles of like 

nature but differing in size only. The growth rate so derived by B & W is 

incorporated with LSW-hydrodynamic model of particle coarsening and the 

asymptotic size distribution determined as a function of the particle volume 

fraction , is in agreement with earlier workers. It is shown that the 

3 

diffusion controlled coarsening the basic kinetic form r is proportional 

to kt, where r is the mean particle radius r , t is time as usual, and k 
av av 

is the rate constant a function of. the volume fraction . 

Bixon and Zwanzig [19a] reexamined the diffusion problem in a median 
with static traps (sinks), and obtained the first order in density. They 
obtained a long time tail solution for its steady state situation included 
more complicated dynamical processes that appear to lead to contributions 
of higher order in the trap density. The work of B & Z was extended by 
Kirkpatrick [20] in which the time dependent transport situation is well 
discussed by adding a class of density corrections which are divergent for 
long times anc that when these divergent terms are resummed they modified 
the loves: order result of E 4 Z. 

1.5. Works by Marqusee-Ross and by Kahiveit 

Marcusee and Ross [ 2 1 a j presented a nev derivation for the las: stage o 
phase separation in the kinetics of a first order phase transition, precipi- 
tation, where Ostwald ripening is the dominant mechanism. They used a time 
scaling [2ia,c,dj technique and derived the power law time dependence and 



distributions for the size of the particles of the new phase. The 
derivation classifies and corrects prior size of' the particles of the new 
phase. Equations are derived for the corection terms to the distribution in 
power laws. The derivation classifies and corrects prior work. In a 
succeeding paper [21b], M & R extend the theory of Ostwald ripening to 
Include the dependence of the volume fraction [ 14-16, 18 ,21b] of the minority 
the precipitate (see Appendix A), the multiscattering [ 18-20,21 b,d] approaches 
(see Appendix C) , and the power laws of the time dependences for the late 
stages of phase separation, Ostwald ripening. 

The work of LS and Wagner (the so called LSW-theory) has been criticized 
by Kahlweit in a series of papers (22a, b]. He claims that their asymptotic 
solutions are not those of the last stage of phase separation but rather those 
of earlier intermediate stage. A 3-dimensional schematic representation 
[22aii] and a 2-dimensional one (22b, c], in distribution of particule number 
density on particule size, demonstrate his idea very clearly. 

Marqusee and Ross [21a] claim that LSV's solution are correct and that 
they are the leading terms in the expansion of long time solution. As for 
the contradiction in fixed total mass of precipitate yet allowing the 
monomer concentration to change, K o t: snow that the assumption does not 
affect the lowest order tern, in the expansion. 

1.6. The main goal of this work 

Tne present paper is intended tc reduce the mathematical intricacy or 
all the pioneers' theories (LSW-theory, MLSW-theory ) in order to simpiiry 
them to an understandable version. Several appendices are added for this 
purpose. Some discussions are also included for stimulating new ideas for 
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Following this introduction, we emphasize the fundamental concepts 
(Section II) and asymptotic distribution for particle sizes (Section III), 
which obeys an equation of continuity in R space (Section III, Eq. 3.1). Key 
points for approximate solution of complicated process, experimental efforts 
following theories, and the theory applicability to alloys are discussed 
(Section IV). A comparison of related works on solution of diffusion equation 
withand or without sink/source term(s) in time period of four decades is 
presented (Section IV, Table 2). We evaluate the growth law (appendix A), 
deduce continuity equation and its Laplace transform (Appendix B) , demonstrate 
multiscattering representation and its Fourier transform (Appendix C) , and 
6how example for solution of Fick's 2nd law (Appendix D). Two figures contain 
in Section I, one table and seven figures are in section III, and one table in 
Section IV. All these effortcs are for simplifying the theory and drawing the 
distance between ideality and reality more closer. 

Most of the related works containing in literature are for systems with 
spherical sinks (precipitate particles). For more complicated system Odijk 
and Lekkerkerxer [47] have recently reported a theoretical approach on the 
phase separation for a solution of bidisperse rodlike macromolecules. It is 
worthwhile to pay attention tc such a system (see Section IV, SuDsecrior. IV. 6). 
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II. 1. Diffusion equations 

Consider a general equation for diffusion, 

dc/dt = D v 2 C + A (2.1) 

o 

in which C is the concentration of a supersaturated solution, D is the 

o 

diffusion coefficient and A is the sink/source term. If replace C by T (the 
temperature), D^ by XV^s (the thermal-conductivity divided by density-times- 
specif ic-heat ) , A by A'/ Ps, then the expression becomes equation of heat 
conduction [-23a]. 

For steady state case,^C/^t = 0, Eq. (2.1) reduces to 

v 2 C = - A/D (2.2) 

o 

Eq. (2.2) is called the Poisson's equation. In case A =0, no sink/source 
and is also in steady state, then, Eq. (2.2) reduces to 

v 2 C = 0 (2.3 > 


which is called the Laplace's equation. Another case is, no sink/source 
bo: is not unde: steady state, from Ec ; . (2.1) one obtain? 




D v~C 
o 


( 


'v 

J 


Eq. (2.?) is called the ecuatior. of continuity. In fluid dynamics, if 
consider v as the velocity of the fluid of density , then V= ^ v is caliso 
the flux density and represents the total flow of fluid per unit cross section 
in unit time. If no fluid is created or disproved vrithin the small volume 
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element d"£, this loss of mass (if any) must equal ~(<t 7«) t)d T , 60 that 


v.V * - df/dt , (2.5) 

Eq. (2.5) is called the equation of continuity l23a] in fluid dynamics. 

If replace p by C and V by j , j=-D Q v C, then Eq. (2.5) becomes Eq. (2.4) 
the equation of continuity. If consider steady state and without sink/source, 
Eq. (2.5) reduces to (2.3), the Laplace's equation. Laplace's equation and 
the equation of continuity are very important in solving diffusion equation 
problems. 

For the last stage of phase separation, when the nucleation is 
neglible, the distribution of particle size obeys a continuity equation in R 
space, where R is the radius of the partice, 

0 n(R,t)/^t + <*>{V[R,C(t)]n(R,t)}/^R « 0 (2.6) 

Where n(R,t) is the number density of the precipitate particles, and 

V[R,C(t)j = dR/dt 

is the growth rate of a precipitate particle^, here C( t ) is the average 
monomer concentration [21a], 


Clt) = C__ - change c-f concentration due to precipitation 

(2.7) 


= C. - j j (4/3) TTR'/v Jn(R,t)dR 
in J < m 


V [ R , C ( t ) j = dR/dt 


= ( Dv / P. ) ! C v t < - C { R) ] ( - • 6 ) 

m eq 


where C. is the initial monomer concentration and v is the molecular 
m m 


volume of the solute. 
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If the solution is not stirred (the diffusion limited case), the first 
step is to transport C to the interface (by diffusion), and the second step 
is to incorporate monomers into the precipitate particle. Therefore the 
quasistationary solution [9b] for diffusion limited growth (dl-case) of a 
spherical particle is 

V dl [R,C(t)] = [C(t)-C (R)]Dv /R (2.9) 

eq m 

For the precipitation of crystals from a stirred liquid solution (the ik- 
case), the growth rate follows a first order law [9b] as 

V lk [R,C(t)] = kv (C(t)-C (R)J (2.10) 

m eq 

Where k is a constant and in dimension cm^molecule *sec *, 

k = D/£ =ifc) (j£) is Lifshitz-Slyozov's notation [9a]), J is density, 

dR/dt = V [ R , C ( t ) ] > 0, if C(t) > C (R) 

eq 

< 0, if C(t) < C (R) 

eq 

Therefore, we may call [C(t)-C „ ( R) ] the driving force for crystal growth 
or dissolution. 


II . 2 . LSW theory and Kahlweit ' s deduction 

According to LSW- theory, the size distribution asymptotically approaches 
a time independent shape, 


n(R,t) 



for t — - oc 


where ^=R/R c , a dimemsionless length, and R is the critical radius. 
This result follows from the assumption that the ageing rate , d^/dt, 
defined in the particular case (22aii) according to LS'W-theor y [9], 
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d§/dt - d(R c )^/dt -*■ b >0, for t ~y *o (see Fig. 1) 

where b=2kv 2<TC(< x5 )/(k_T) , in which k_ i6 the Boltsmann constant, 
m Is B 

the other constants have the usual meaning. Kahlweit's derivation (22aii) 
follows from the conclusion that the ageing rate defined in this particular 
case by d^/dt = d(R ;/dt. Since according to LS-theory [9a, Eq.2], 

dR/dt = O0/R)(A - *Vr), ^ = ^/R , <*/R = C(R)-C(^°), so that 
dl/dt = 2R dR/dt 

= 2£(^- oC /R) 

= 2/C[(oc/R c ) - C(R) + C( oo ) ] 

Since as t —yoo t C(R) —^C(°"), therefore 

di/dt -> 2ko6/ R = 2k[ 2 CT v C(^)/kT]/R = b > 0 

c m B e 

in which at? = D/j has been replaced by k, the rate constant for this particular 
case, precipitation of crystals from a stirred liquid solution [22aii]. 
However, Kahlweit (22ai) showed that the latter conclusion is incorrect. 
Instead, d 5/dt may be seen as a measure of -dC/dt as shown in Figure 3, 
increases rather rapidly after the nucleation period and then passes through 
a maximum curing the early stages of precipitation to slowly approach zero 
for C= C*. 

II. 2. Sind den; 1 1 y and competition of sinks in preclpl t at ion 

Consider the total number or molecules (or moles) in one sink change; 
vi. z r. £ i Qc S.S 

d( 4 7rF."/3v )/Jt = k(R) P 
m ( 

« 

= k(R)[C av -C eq (R)j 
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4VR 2 dR/dt 


(see Appendix A) 


The above expression is for one sink only. In competition with other sinks 
the expression [21b] will be 

4F R 2 dR/dt = 4*D o R[ l+Rj4T^<R>][C av -C (R)] 

dR/dt = (D /R)[ 1 +rJ4^7<R>] [C -C (R)] (2.11) 

o * \ av eq 

II. 4v Particle size plays an Important role in precipitation competition 

The driving force for ageing of the precipitate is originated from the 
surface free energy difference between different sizes of molecules. The 
surface free energy is related to the size of the droplet for a liquid-vapor 
system in the following way: The vapor pressure of smaller droplet will be 

higher than that of the bigger droplet as follows, 

In (p 2 / Pl ) = (2"/ M/f d RT)(l/a 2 - 1/a^ 

= (2 i v /kT) ( 1 /a - 1/a,) (2. 12) 

m 2 1 

where ^ is the density of the condensed phase precipitate, "/ is the 

surface tension, p^, p^ are the vapor pressures of the droplets with 

radii a, , a respectively. H/RT is the molecular weight divided by 
i 2 . 

RT, gas constant times the absolute temperature. This fact or principle 
favors the growth of larger particles and derives the precipitate to decrease 
its surface/volume ratio, that is growing bigger in expense of dissolving the 
smaller particles. 

Similar deduction can be applied to solid-liquid phase transitions. A 
number of observations of the growth of solid particles in a liquid medium at 
the expense of smaller particles have previously been reported [7c,d], 
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Greenwood [7c] proposed a general analysis of particle growth using 
Thomson-Freundlich equation [7f]: 

ln(S/S a ) - 2M<f7(RTf d a) (2.13) 

Where S is the true solubility, S the solubility in contact with a particle 

d 

of radius a, (T the interfacial tension between solid phase and liquid phase, 
and M/RT has the usual meaning. Bigelow and Trimble [7d] repeated Ostwald's 
experiment [7e] on sulfur and perhaps led the name Ostwald ageing [9b] and or 
Ostwald ripening [18, 21a, b] for the last stage of precipitation in binary or 
multicomponent systems. 

For a special particle and an ideal solution, the R depencence of 

C (R) is 
eq 

C (R) - C (*») exp(oC/R) 
eq eq 

= C (^)(1 + <* /R) (2.14) 

eq 

where = 2^ \ /k T [21a, b], 
n> 15 

Conventionally, a set of reduced variables are used as follows: 

r=R/oC , T= t/T , T G ~=od i ’/Dv C (** ) and 

m eq 

T ik=oc /'[kC CF \ T >=; C(t)-C (**=•))/€ (2.15) 

eq ‘ ec eq 

in which T=(R )' > /o t*b in LS-notatior. [9a]. 

co 

Briefly, there excsts a critical radius P. which is the particie 

c 

radius in equilibrium with the surrounding solution, and R is the effective 
particle radius [9b] defined as 

R = 2 dV/dO = 2 d[(4/3)7rR 3 ]/d(47R 2 ) (2.16) 


quality 
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The statistical theory of number density as a function of particle size 
first established by Smoluchovski (24c] as a linear function of time t iS 

1/n = l/n Q + const. t (2.17) 

where const = 4TOA, in which A is the radius of action, D is the 

diffusion coefficient, and n and n are the number densities at time 

o 

t and t respectively. Smoluchovski also reached a mathematical 
deduction formula for the number of monomers in cluster in such a way 

= n o (t/T) k_1 /(l+t/T) k+1 (2.18) 

where n^ is the number density of clusters, in which each cluster contains 
k monomers, and T is a constant which equals (4TTDAn o ) *. 

11.5. Mass transfer in crystal growth 

The interface of a crystal grain represents not only a sink for 
matter but also a source of both heat and matter. The rate with which 
crystal building blocks are incorporated into the solid is governed by the 
difference of their chemical potential in the two phases in contact. 

Though for the time dependent position and shape of an interface one must 
simultaneously consider the conservation of mass, momentum and energy for 
the system ir. question, however, one mav simnlifv the problem by reasonable 
assumptions (such as steady state approximation, isothermal condition, etc.). 

A critical discussion on diffusive mass transfer in crystal growth was 
given by Westphai and Rosenberger [27]. The discrepancies between predicted 
and measured growth rates were attributed to the formal description which 
ignored the limited analogy between heat transfer and mass transfer and the 
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Nernst -Brunner unstirred film model. However, the intricacy of mass transfer 
in crystal growth can be simplified by proper approximations. Drastic 
simplications can be made for the crystal growth continuity equations - rate 
of gain in mass, in momentum, in internal and kinetic energy [27b, Eqs. 5.1 - 
5.3). As for convection may be neglected when the crystal grains are smaller 
than about lO/i, because the crystals move with the fluid [28]. Under such 
circumstances one would expect the growth rate to be insensitive to stirring. 
Furthermore, the diffusion fields of the particles have been found to be 
independent of one another if the average distance between them is larger 
than about 10 times their diameter. Under such circumstances diffusion 
controlled growth has been observed [29,30]. , 

Phase equilibria and chemical equilibria depend on temperature, thus 
a temperature difference can cause a solid to dissolve in a fluid at one 
temperature and deposit elsewhere at a different temperature. Gradients of 
many types cause relative motion of the components of a mixture with 
respect to one another. Hurle et. al. pointed out [31] that transport takes 
place whenever a difference in chemical potential exists, such as between 
stable and metastable phases of the same material. 

By experimental measurements one may reduce the intricacy of the 
diffusion mass transport problem by imposing an isothermal condition. Aisc 
one may focus attentions to the reference frames chosen for the flows of the 
entire system. 

For studies or. multicomponent isothermal diffusion it is customary tc 
describe the transport matter by means of flow equations which are extensions 
of Fick's first law. These equations express the flow of current J.. of 
component i as a sum of concentration gradients each multiplied by a diffusio 
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coefficient (named as practical diffusion coefficients) of one type. In 
theoretical work the corresponding phenomenological equations are usually 
written as sums of forces, i.e., the negative gradients of chemical potentials 
each multiplied by a diffusion coefficient D'^ of another type [32]. The two 
sets of diffusion coefficients dependent on the frame of reference selected. 
Onsager's reciprocal relations [33] hold between certain of the D' ^ and the 
testing of these relations by using data for the together with certain 
thermodynamic data. Correct tests of the Onsager reciprocal relations 
require proper specification of the reference frame. Usually, the solvent 
fixed frame and the volume fixed frame are used. The former is used when 
the chemical potential gradient of the solvent is eliminated from the flow 
equations, and which is used for computing D' ^ from , and then for 
evaluating these from those for the volume fixed frame. The volume 
fixed frame is important to the experimental measurements of diffusion. 

For the solvent fixed or mass fixed frames some auxiliary data such as 
partial specific volumes are required for determination of partial diffusion 
coefficients. Kirkwood [32] considered, for certain frames of reference, 
as is the problem of measuring practical diffusion coefficients correctly 
where the partial specific volumes are dependent on concentration and the 
solutions therefore exhibit a change of volume on mixing [32]. 


In addition to isothermal condition, the phenomenological description 
of diffusion in isotropic mixtures has been developed from the expression 
for the entropy production and the relation between the diffusion 
coefficients in accordance with Onsager reciprocal relations [ 3 4 j . 


11.6. Entropy production and Onsager reciprocal relations 

The phenomenological description of diffusion in isotropic mixtures is 
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derived from the expression for the entropy production and the relation 
between diffusion coefficients, resulting from the Onsager relations [34a]. 
The entropy production of the vectorial irreversible phenomena can be 
written in the form 


O' - Z j k . X. 

k=l 


(2.19) 


where are the components of fluxes and are the forces. In 
isotropic media 


J k • , k-1 

j-i 3 3 


= 1 , 2 , 


( 2 . 20 ) 


In the absence of magnetic field, the Onsager relations [4b, d] state that 


j = Lj k » j , k— 1 ,2 n 


( 2 . 21 ) 


The system described by Eq. (2.15) can be solved for the X ^ , the new scheme 
of phenomenological coefficients again being symmetrical, i.e., 


X, = X C 1 ^, 


k-1 


jk k’ 


, "I _ T "I 

i — - L. . 

J K kj 


( 2.22 


The entropy flow J ant entropy production (j car. be deduced from the 


thermodynamics equation Td5. = dE t 


PdV - X - A3 i dc - • lr d 


t demonstrates that 


the entropy prodection contains four different contributions [35, Eq. 21 


1) heat conduction (due to temperature gradients). 
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2) diffusion (due to concentration gradients), 

3) viscous flow (due to velocity gradients), and 

4) chemical reaction. 

For systems without heat conduction and without chemical reaction the 
contributions reduce to two. 

II. 7 . Thermodynamical argument on crystal growth 

Since by controlling the heat removal from the fused metal the 
metallugist can give the casting the requisite structures, Ivantsov studied 
the growth kinetics of different crystal structures [36a], When the crystal 
is small and the growth rate very small, the system will be almost 
isothermal, i.e., the effect of uneven temperature distribution will be 
vanishingly small. 

In addition to the diffusion process, the molecular kinetic process in 
crystallization under supersaturation is equally important [36b], especially 
for controlling the crystai-grain sizes [36c], The fundamental aspects of 
crystal growth concern: a) the atomic processes of phase change, b) surface 
energy, and c) heat and matter transfer, a) and b) influence the geometry o 
the interface, c) determines the temperature and compositions of crystal 
growth [36d]. The interface roughness is a measure of entropy of fusion 
6ince the interface being rough for metals and smooth for most nonmetais, an 
the entropy of phase change is small for rough surfaces and large for smooth 
surfaces. According to this fact one may find, for some reason, compounds 
which may have low entropy of fusion freeze like metals [36e]. 

The spontaneous precipitation is a process toward a state of low 
entropy. However, from the point of view of statistical mechanics, all 
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molecular processes must be reversible. Besides, there i6 a theorem in 
dynamics due to Poincare, which states that in a system of material particles 
under the influence of forces which depend only on the spatial coordinates, 
must, in general, be able to return to the initial state to any desired degree 
of accuracy [24bii]. 

In accordance with Boltsmann's point of view bssed on probability, 
Smoluchovski's theory of fluctuation in molecular concentrations [24] allows 
us to bridge the gap between the region of the macroscopically irreversible 
diffusion and the microscopically reversible fluctuation. 

thermodynamically a spontaneous process can occur when the free energy 
change is less than zero even if the entropy of the system is decreased. The 
Gibb's thermodynamic potential, i.e., the Gibb's free energy, can be deduced 
for the precipitation in supersaturated solution as follows: 

1) dG = -Sdt + VdP + YriA + ^?G/^n.dn. 

i 1 1 

dG n = -SdT + YdA + J G.dn., (2.2^) 

P ~ i l 

during the latent perioc, dT > 0, so -SdKO, where S is the entropy 
of the entire system. The 3rd term of Eq. (2.24) is < 0, since 
dn^ < 0 for the solute while precipitation occurs anc n of 
the solvent is unchangec. nance dG<0. Note that though the 2nd term 
i dr. should be >0 (the interfacial tension Y>0 , dA >0) for one crystal 
growing, however for the crystal growth in expenditure of scalier 
particles a summation of y. dA. should be used. As long as {dA..} are 

negative except the one of the crystal which is growing, therefore the 
3rd term has no significant contribution to the free energy change. 
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Consider the phase change being isothermal. 


dG - d(H - TS) + i dA + 2. G; dn^ 

i 

« dH - TdS + Y dA + £ (F; dn c (2.25) 

i 

in which dH < 0, -TdS > 0 (sine dS<0 in precipitation), the 4th term 
is less than zero, and the 3rd term has small significance as discussed 
above, therefore we expect dG<0. 



■-StoStS 


III. ASYMPTOTIC DISTRIBUTION FOR PARTICLE SIZES 


U.II.1. Fundamental equations and reduced variables 

Define n(R,t) as the number density of particles with radius R per unit 
volume of solution. N(R,t) can only change by growth from nuclei or 
-dissolution of precipitate particles. It must obey a continuity equation in 
R space in the specific dynamic system, 


^n(R,t)/^t + ^/<3 r{V[R,C (t)]n(R,t)} = 0 

3 V 


(3.1) 


Where {V[ R.C^C t) ] n(R , t ) } is called the flux density and represents the 
total flow of fluid per unit cross section in unit time [23a]. It is also 
assumed that the particles are stationary and nucleation is neglected, i.e. 
the 6ink/ source term in Eq. (A. 1) of Appendix A is nearly zero and is omitted. 
The precipitate-particle density at time t is N(t) (see Fig. 3), which is 
conventionally termed (t), and expressed as 


N(t) = f n(R, i 
'0 


,t)d?. (3.2) 

une can see from Fig. 3 

•o 

n(R,t,)dR = Nvi.) and i n(E,:JdR = N(t ) 

0 1 ; '£ A 1 

Tne (n,t) plane, it Fig. 3, illustrates the course of development of 
K— U 

tne mean concentration of monomers with time [22aiij, which is related to 
cavity theory [36,45.' and theory of fluctuations [24bii]. The (R,t) 

D“v.- 

niane shows the trajectoory of the projections of critical radius of the 
clusters changes with time. Tne planes parallel to the (n,R) ti .Q plane 
represent the size distribution as it changes with increasing time. 
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The distribution function can change by: 1) growth of particles, 

2) dissolution of particles, and 3) the creation of new particles due to 
nucleation. For last stage of phase separation, in case the nucleation rate 
is negligible, the distribution obeys an equation of continuity in R space 
as shown in Eq. (3.1) above. 


As mentioned in Section II, by introducing reduced variables "r" or "a", 
T , and cT (T) [21a, b], the continuity equation in reduced variables [21e] 
becomes: 


•^TF(a.T) +^/<)a{F(a,T)V[a,<f(T)]} - 0 


(3.3) 


with the growth law [ 2 1 f ] as: 

V[a,<rTO] *= a X [<T(T) - 1/a], A=1 for dl-cases 

=0 for ik-cases (3.4) 

The dimensionless reduced variables are defined [ 2 1 g ] as: 


a=R loC , T =t/T, O' (T ) = [C(t)-C (*>)]/C («») 

eq eq 


(3.5) 


III. 2. Decay function in supers a tur at lor. 


CT ( X ) is the supersaturation at t >/ 0 , at t = 0 (T(0) is called (J'T . 

f in 

As time goes or, as Ostvald ripening begins, ) decreases, so that it is 


a measure of the decay of the supersaturation, i.t., 


O'CT ) = (T -X } a J F(a,T)da 
m 

The constants in Ecs. (3.5) and (3.6) [ 2 i h ] are: 


(3.6) 


oL = 2 V v /k I , 

m c 

T ik = 06 /kc («), T dl =o 6 2 /D v C (*>), and 
eq o m eq 
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(3.7) 


K- - 4T^ 3 /l3v B C e{j (»)] 

where 7 is the surface tension and C (P°) is the concentration in 

eq 

equilibrium with a macroscopic particle. 

The ageing rate and the time development of the size distribution of a 
precipitate have been analyzed both experimentally and theoretically [22] on 
micelle association and dissociation in surfactant solutions. Kinetic 
equations have been given a form which suggests an analogy with heat 
conduction. Experimental methods, mostly in relaxation times, have been 
deduced for the net distribution of the rearrangement among the micelles 
[37]. 


Equations (3.3) to (3.7) are the starting point of the analysis of the 
last stage of phase separation [21a). In order to find the asymptotic solution of 
the partial differential equation, the continuity equation, (3.3), one has 
to analyze N(t) anda'(T). After the rapid growth of the newly formed 
nuclei, the number of the precipitate particles N(t) and the supersaturation 
O' (T) are known to decay slowly. The number of colloidal particles decreases 
slowly. The supersaturation (f(T) decreases slowly from (T<1 to approaching 
zero. The total mass of precipitate increases slowly. The average radius 
<R> increases slowly. All these imply that no exponential law of decay but 
power laws [ 2i a . b. . 

III. 3. Power expansion and scaled distribution function 

Since ^F/r't =?=■ 0 (ir. the sense of fluid dynamics, it is not 
incompressible) , therefore 

v.V-=^=0, V = V[a,<T(T) ]F(a,T) , i.e., 
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v.V “^/3a {[V(a,<T(T)]F(a,f)} , for a spherical-function-aolution. 

A power expansion which is sufficient for the determination of the properties 

* * 

of spherical function defined in the interval 0 < a < a , where a 
is a characteristic constant, which can be used as the asymptotic solution 
of the differential equation. In such cases the differential equation can be 
integrated "asymptotically" for very great values of the variable [23b], 
so that the solution can be interpreted asymptotically. 


In extracting the lowest decay from the distribution function and write 
it in the long-time limit [21b], Marqusee-Ross's scaled form for the 
distribution takes the following form: 

F(a,T) =T' y F o (z); z=aT~ X 

or, takes the form as in their preceding paper [21a], but setting r = a, 
one obtains: 

, * 

F(a,r) = 7 ZT (z) (3.8) 

n-0 n 

z = aT X and <3 /c> a = T X <? id z (3.9) 


Substituting Eqs. (3.8) and (3.9) into (3.3) and (3.4), 
continuity equation and deduce the growth law in reduced 


showing in (3.10) and (3-11) respectively (see below). 
^ F(a,Z )/3?c + sfea [F(a.T) V[a, r(X)U = 0 


V[8,CT(T) j 



(T ~C-n 


1 /z)7 X j 


one may solve the 
variable forms, as 


] 0 ) 


v 


II) 


A* 

In which G~ ^ T '' r. has been introduced. The distribution of 

, n 

n= 1 


long time limit in its scaled rorm is: 
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F(a, T) * *r‘ y F 0 (z), 

As in the method of series expansaion, (3.10) can be solved by 
substituting (3.8) and (3.11) into (3.10) with the aid of (3.9). By 
collecting the coefficients on both sides of the ultimate equation, for 
a certain n, the following rules hold: 

a of F (z) on LHS = a' of F (z) on RHS 
n n n n 

b of dF (z)/dz on LHS = b' of dF (z)/dz on RHS 
n n n n 

where a , a' are coefficients of F (z) on LHS and RHS 
n n n 

respectively; b , b' are coefficients of dF (z)/dz on LHS 
r n n n 

and RHS respectively. 

Using these rules one finds that (2 +X)x -1=0, therefore 
x = 1/(2 +X) 

= 1/2 for X = 0 ik-case 

= 1/3 for A= 1 dl-case 

Consider two values of n, n=G and n=i , one obtains: 

4x F (z) + xz d/dz F (z) = d/dz 't \(T.- l/z)F (z) j (3.12) 

o o j o 

3x F ( z ) + xz d/dz F . ( z ) = d/dz { z ) F . ( z ) 

^ i a 1 

-JZF(z)} (3.13) 

Z o 

The operator d/dz, on the RHS in both Ec . (3.12) and Eq . (3.13), may 
generate dt^(z)/dz and F^(z) terms which car. be combined with the LHS proper 
terms. In this way one obtains the scaled continuity equations, and scaled 
growth rates £)(z, Sf ) as follows. 
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-F (z) + d/dz {U)(z, O'.) F(z )} - 0 (3.14) 

O 10 

-F^z) + d/dz {(3/4W( z , Op F x (z)} 

+ d/dz {^>(z, 0^) F q (z) } = 0 (3.13) 

U>(z, (Jp « z _X [ - l/z]/3x - z/3 (3.16) 

*Hz, Op - z _> (T/4x (3.17) 


Here dJ(z, *T ) are the "scaled growth rate" expressions, 
n 

III. 4. Scaled growth rate as function of scaled radius 

Plottings for cJ(z, Tp versus z with different -^-and ^ are shown 
in Fig. 4 and Fig. 5. The calculated data for the plottings are listed 
in Table 1. 

Figs. 4 and 5 show that as u)(z, <7p = 0 and d<d(z, cp/ dz = 0, one 
obtains the optimum values of z, as z q , and 0^ the decacv function of the 
supersaturation, for different values of X and x, i.e. , 

0, x = 1/2, z = 2 1/2 = cr 
o 1 

-X = l, X = i/3, z = ( 3 / 2 ) 1 7 3 , = (S/4) 1/3 = z 2 (3.18) 

0 J o 

The general font of tne distribution function, F (z), in the last state 

n 

of phase separation, can be solved by integration using separation of 
variable and other techniques. For instance, consider Eg. (3.14), by 
separation of variables, one obtains: 

dF (z)/F (z) + d uKz,<r)/^( z ,<r ) 
o o 11 

= [<^(z, ^)) -1 dz (3.19) 
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(3.19) 


which is d ln{F (z) tJ(z f O')} « («>(z, T!)] * dz 
o 1 1 

III. 5. Scaled Distribution in Zero Volume Fraction of Sinks 

After integration Eq. (3.19)' the scaled distribution function, F q (z), 
is obtained as: 

F q (z) «= const[ A)(z, <^)] 1 

( Z -l 

x exp J (zJ( z' , CTp ] dz', z <= z q ; 

=0, z > z q (3.20) 

Substituting (3.16) into (3.20) for -X=0 and A=1 respectively, and using 
the proper values of z^ and 0^ (i.e. using Eq. (3.18)), then by tedious 
manipulations, one obtains: 

1. For \= 0, x = 1/2, z = 2 1/2 , CT = z = 2 l/2 

o 1 o 

l^Cz, 1 = -3z/(2 i//2 - z) 

F (z) = (C CT /K ) z (2 i/2 -z)" 5 
o o in 

• - / / _ 1 / 2 .. , - 1/2 
x exp i — _>z ■ (2 -z)j, z <= 2 ; 

= 0 , z > 2 * ; ‘ (3.21) 

2. For K= 1, x = 1/3, z q = (3/2) 1/j , CT = (9 /4) 1/3 * 

. - 1 7 . 'i 7 

!iJ(z, J, ) j ' = -3z“ / 1 z' - 3z ~z + 3j 

I c 

F (z) = (C CT // c)z 2 (2z +z) 2y/_ (z -z) 
o o m o o 

x exp! — z / ( z o “z)j > z <= (3/2) 1 ' /3 ; 
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2 > (3/2) 


(3.22) 


0. 


1/3 


The above results are obtained by using various integrating techniques, 
including change of variable and rational fractions [23c]. Eq. (3.22) is for 
the dl-case and with the assumption of negligible nucleation, i.e., =0. 


By use of the imposed normalization condition for F q (z) [21a], 


or 

in 


z F (z)dz 
o 


(3.23) 


one may substitute F q (z) for z <= z q , in expression (3.22), into (3.23), to 

obtain the constant c . Numerical calculation method may be used for 

o 

evaluating C q . Since J ust the normalization constant for 

3 

z F q (z). The justification may be made by rearranging Eq. (3.23) 
to the following expression: 


or uc , p 


z F (z)dz 


4 


= const / z f (z)dz 

J o ° 


= const/ c 


(3.23)' 


in which c 1 = ; z J f (z)dz has been imposed, so that 

O -V o 


( 07 i^- )c = const 

That is to sav, if the integral car. be evaluated, c * is obtained, 

o 


so is the constant ir. Eq. (3.23)'. 


For expression (2.21) of z <= z , evaluation of c may be evaluated 

o o' 

by change of variable several times, then the integral may be reduced to a 
standard form named "exponential integral", so that proper mathematical 
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table can be used [23d] 
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III. 6. Scaled distribution in nonzero volume fraction of sinks 

With <^=^0, correction term should be added, it can be written [21b], 
after another tedious manipulation, for instance, for A- = 1 , dl-case, 


£*>(z, (Tj) * (l/z)[ O' - l/z][l + z6 ] - z/3 
F_(z) * C_z 2 (z_-z)”°^(z+3/z_ 2 ) P 

.1/3 


o o 


x exp[-Jz/(z o ~z)] , z <= (3/2) 


0, z > (3/2) 


1/3 


(3.24) 


Here OC = 


2 + (3z ^ + 18z S/(z 3 + 3)' 
o o o 


e 

s 


1 + 27/(3 + z 3 ) 2 
o 


3z ^/(3 + z 3 ) 
o o 


(3.25) 


and 


6 - 


1/2 


rtC 


f 

4U = (4TT /3) ! R 3 n(R, )dR 
Jc\ 


(3.26) 


Here is the equilibrium volume fraction of precipitate in solution (or 

in melt) after given quench, or eouivaientiv = ;C. - C ( p ")]v [21b] . 

' ir. eq m 

The effect of equilibrium volume fraction on the scaled distribution 

function can be seen from a set of plottines. F (z) vs z in di-case for 

-o' 

different f ec (f eo =0, 0.005, 0.015 , 0.035, and 0.075) as shown in Fig. 6. 


It is interesting to note that. 
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R [4ir^<R>] 


1/2 


1 / 2 

a [4TT^<a>] , when replace R by oL a, and ^ by £ /oC 


T, .3 


1 / 2 

z I3 4eu a 0 ^,3 » w h en the following conditions are imposed: 

1) n(R,t)dR = [3<|W(4'ncd 3 )]F(a,r )da. 


I 


oo 


2) ) z F (z)dz = 1, and 

'0 


3) a = z“£ x , <a> = a Q T X , ^ (a) = T *• 

By numerical calculation method, plottings for z q vs t~ -T< reduced 
coarsening rate) vs <b , and a , (J" , P vs & are shown in Figs. 7, 8, 

IPO O llo > 

9 respectively. From Fig. 9 one can see that, as slowly increases with 
time, the average radius a^C^*, ) increases slowly, the supersaturat ion 
& j ( 4 so ) decreases slowly, but the density of particles or droplets of the 
minority phase ^(^po) decreases faster than the change of degree of 
supersaturation. It can be deduced in this way, even if the supersaturation 
reaches a minimum ^ still changes fast because the expenditure of 
smaller particles in growing the bigger ones should drastically reduce the 
number density Therefore the following power laws may be considered 

as the time rate as if instantaneously being kept constant. 


- 1 


f- 1 "- 


anc 
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Table 1. 

u)(z. 

-X 

CTj) - z ( 

C^-z *)/3x - z/3, x*(2+X.) * 


z 

x=l/2,lk-case,fc<z, <jp 

X«l, x*l/3, dl-case, 

l 

4J(Z, Op 




-2 1/2 

-2 

Cj-1 “(9/4) 1/3 

-2 1/2 

*2 

1/2 

-0.83 

-0.56 

-0.17 

-2.17 -2.05 

-1.34 

-0.17 

1 

-0.33 

-0.06 

0.33 

-0.33 -0.02 

0.08 

-0.67 

1 / 3 

(3/2; | 




0 



2 1 / 2 


0 





3/2 

-0.28 

-0.002 

0.39 

-0.28 -0.07 

-0.002 

-0.40 

2 

-0.33 

-0.06 

0.33 

-0.42 -0.26 

-0.21 

0.08 

3 

-0.56 

-0.28 

0.11 

-0.78 -0.68 

—0. 64 

-0.44 

4 


-0.56 

-0.17 

-1.15 -1.07 

-1.04 

-0.90 

graphs, 

(See Fig. 4 , 4)= 

(See Fig. 5, Id 

z=0, as 


60 vs z 

as z ■= 
o 

cr.2 1 ' 2 

J i 

.) 

z =( 3/2) 1/3 , < = (9/4) 1/3 .) 
0 I 
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IV. DISCUSSION AND CONCLUDING-REMARKS 


IV . 1 . Concentration gradient as driving force in diffusion 

Consider one of the precipitate particles (sinks), around it there is a 
concentration gradient for the other species, and that the rate of flow 
of solute molecules in this concentration gradient is governed by Fick's 
law of diffusion. The Smoluchowski' s theory on the kinetics of colloid 
coagulation as a diffusion-controlled process using Fick's law of diffusion 
[24a, 38] led series of investigations and applications in this field both 
theoretically and experimentally [39]. 

The mathematical solution for the fundamental differential equation 
governed by Fick's 2nd law drew attentions of a great number of physicists, 
chemists, and scientists of related fields [39,40]. Here we deal with the 
situation around one sink particle in very law sink concentration. In this 
case, the concentration gradient for the unprecipitated solute in solution 
is set up and the rate of flow of monomers of the solute into (or out from) 
the sink (or source) can be solved by imposing proper boundary conditions. 

For instance, in case the system is increasing supersaturation, i.e., the sink 
shrinks and is becoming a source (giving up monomers insteac of absorbing 
monomers), tne concentration gradient of the monomer is negative, vC(r,t) d 0. 
As the Ostwaic ripening begins, the concentration gradient will be reversed, 
i.e., vC >0. sc that the monomers diffuse into the precipitate (the sink,-. 

IV. 2. Solut ior. of F 1 c k ' s law, com par isor. between selected works 

— > 

For vC > 0, C(R,t) = 0 is a proper boundary condition for Ostwald 
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ripening process both in physical aeaning and mathematical logic. Since in 
Ostvald ripening the monomer whenever strikes the surface of the sink, shall 
be incorporated with the sink encountered. However, if not every monomer 
reaching the reaction radius incoporates with the sink, an obvious 
modification is to be assumed (38). 

In Table 2, we outline several representative works in this field, in 

different time period (1949 to 1984), from which techniques in solving the 

2 

fundamental partial differential equation, ^C/^t = Dv C + A, with A = 0 and 

or A=^0, in accordance with imposed boundary condition (B. C.) and initial 

condition (I. C. ) , are worth learning. For instance, the boundary and 

initial conditions for the last stage of precipitation (Ostwald ripening) 

are C(R,t)=0 and C(r,0)=C o R<r<=*'° respectively in most of the methods 

cited in Table 2. However, there is a singularity at t=0 for the solution 

chosen by Collins & Kimball [38], since they used the technique in solving 

second order partial differential equation on linear flow of heat [41], 

Also, since C(R)=0 can not be correct if not every monomer which reaches the 

reaction radius reacts. Their modification for the boundary condition is 

C(R) = y(< ? C/V r) . We can see that the modified B. C. still holds if 
K 

C(r) = constant so that (cC/<?r) = 0. 

K K 

IV. 3. Appr oxlmate solution for complicated process , example 

Reiss and La Her [42] proposed a very interesting idea: moving 
boundaries, connected with the growth of coilidal particles. For this 
complicated process, the differential equation is difficult tc treat 
mathematically. In considering the boundary being not stationary, the 
boundary value problem is not soluble by means of the familiar 
methods. So that approximate solution techniques were used [42a] . 
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Interference of growing spherical precipitate particles were treated 
by Wert and Zener [43a] using approximate method in a relatively simple 
manner. A transformation law to fit the experimental data was developed 
by Wert [43b], and results of derivations of Zener [43c] were used to 
calculate the mean distance between nucclei [43b], Following are the 
concepts and equations they developed and applied for experimental fitting: 

1. The amount of precipitate formed at any time is proportional to 

C - C(t). 
o 

2. The rate of transformation of the particles of a second phase 
growing in a matrix without mutual interference may be written as, 

dW(t)/dt *= f(t) (4.1) 

where W(t) Is the fraction of transformed material and f(t) is some 
function of time. However, if interference between precipitate 
particles are considered, the rate will be lessened and the decrease 
will be proportional to the fraction already transformed, so that 
the modified equaztion becomes, 

dV(t)/ct = [1 - W(t)]f(t) (4.2) 

3. Eg. (n.2) can be integrated directly and the ir.tegratec font is 

t 

W(t) = 1 - exp- \ f ( t ' ) d t ' ] (4.3) 

' C> 

An explicit form used by Wert [43b] for Eq. (i.3) is 
C(t)/C o = 1 - exp[(- (t/lT) 0 ] (4.4) 
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where ^C(t) is the amount of transformed material, A C( t)“C o ~C( t) , 

C is the initial concentration of solute, T is a time constant 
o 

depending on C q , T (the temperature of the system), and other 
undertermined factors. 

X can be determined by fitting experimental data to log-ln-form 
of Eq. (4. A), i.e.. 


log ln[C o /^l C(t) ] = n logt - n log (4.4)' 

by plotting ln[C Q /^C(t)3 versus t in log-scale, for a certain 
temperature T, the slope is n and the intercept is X n , so that 

T . -1/slope 

= intercept 

Data for the precipitation of Fe^C in ibC-iron at various 
temperatures were used for fitting [43b]. By taking average 
of the slopes for different temperatures they obtained 

n-1.45 


Which can be approximated as n=3/2 as being used conventionally 
[43dj. The slope, n, varies between 1.2 to 1.7 from 43°C to 312°C, 
while X varies from about 12000 minutes to 0.5 minute [43b] . 

The radius of the spherical precipitate particle, S, varies with 
time [43cj according to, 

S = CxC-x (Dt) 17 - (4.3) 

where D is the atomic diffusion coefficient, is the 

dimensionless growth coefficient, in which X is the number of 
dimensions. For shperical particles =5 of course. OC, takes 
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the form: 


oC^ - n J )/(n o - n^)] 1 ^ 2 (4.6) 

Where n^ » concentration of solute in precipitate, 

** concentration of solute in the matrix which is in 
equilibrium with the precipitate, and 

n^= concentration of solute in the matrix far away 
from the precipitate. 

IV . 4 . Experimental efforts following Theories 

Light scattering studies of phase separation in fluid mixtures [lOe-lOh] 
have provided a wealth of information about the complex, highly nonlinear 
processes by which a new phase forms in a system that was originally 
homogeneous. By light scattering and by direct microscopic observation 
Chou and Goldburg [ 1 0 g } confirmed the early conjecture of Cahn and Moldover 
(Ref. 1 of { 1 Og j ) that hydrodynamic effects are important in the late stage 
of the process, 6ince droplet growth is speeded up by fluid flow which in 
turn is driven by the curvature of irregular surfaces which interpenetrate 
the nucleating fluid. Wong and Knobler [10fj, using an ingenious pressure- 
jump technique to quencr. the system, mace extensive light scattering 
measurements of pease separation in isobutyric acid and water. They found 
that when the critical mixture is quenched, the two-phase region scatters 
light in the form of a ring which brightens, sharpens, and decreases in 
diameter as the system evolves toward its final equilibrium state. 

In accordance with Siggia's theoretical hydrodynamic arguments [10a], 
Wong and Knobler [lOh] experimentally found that the crossover wave number 
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corresponding to the maximum in the scattering, k , decreases initially as 

m 

-1/3 -1 

t and at long times as t . Their earlier work [ 1 Of ] shows an 
2 

initial t growth in I(k m ), the intensity of scattering, followed by 

1/2 

a slower growth t , correspondingly Siggia [10a] proposed a crossover 
1/3 

from r _ t . 


THe existence of crossover can be related to the volume fraction of the 
new phase, <p . Wong and Knobler [ 1 Oh ] found that in dilute mixtures 

k m =At“ l/3 with A=A(c f>) for 0.02 <= f <= 0.10 (4.7) 

and they found their case for <f> < 0.02 is equivalent to Lif shitz-Slyozovc 
growth in which A is supposed to be independent of <f > . Wong and Knobler 
[lOh] aalso compared their results with the quantitative nonlinear theory by 
Langer, Bar-on, and Miller [10b] usually called LBM-theory, and with 
computer simulation on model systems, they found close agreement between 
experiment, theory, and computer simulation [ 1 Oh ] . 


IV. 5. Application of theory to alloys 

It is well known, in fluid dynamics, the diffusion constant of a 
spherical drop of one fluid in another of similar viscosity [10a, 26b] is 

D = k.T/5 nrpys (4.8) 

D \ 

Where a is the radius, fi is density, and ^ is viscositv. Consider drops 

t 

of radius R. , diffusion constant D,. , and number density f* . , for 

X -i- 1 

i=i,2, one finds that, according to collision theory, the number of coliosio 
between the two populations per unit time per unit volume [10a] is: 

X 1 2 = 4U(R 1 +R 2 )(D 1 +D 2 )f ] (4.9) 
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If each collision causes one new droplet formation, then one may expect a 
population of drops agglutination at a rate given by (4.9) approches a 
stationary distribution when scaled by average values R, D, and i.e., 
R 1 +R 2 »2R, Dj+D 2 =2D, and with the assumption Cj ^2 “ ^ * then 

tf/H « -165TD(R)R^ 2 (4.10) 

with (4tr/3)R 3 ^>= <f> (4.11) 

Where D(R) is the diffusion coefficient for a particle of radius R and <f>±s 
the volume fraction of particles. 

By using the definition of kinematic viscosity [26c], and 

substituting \l ^ for ^ into (4.8), also replace "a" in (4.8) by R just for 
convenience, one obtains 


D(R) = LT/5 TT R 

D 


(4.12) 


Substituting (4.12) into (4.10) with the aid of (4.11), after integration^ one 
obtains the power laws 


n -1 , „ » 1/3 1/3 

P _ t and R _ <£> t 


(4.13) 


However, for the case of a binary alloy where we have a solid parricle in a 
solid matrix or in very high viscous fluid, the diffusion coefficient at low 
temperatures is estimated [i6t,2ibj to depend on the radius as 


D(R) - 1/R 


(4.14) 


and at high temperatures as 


D(R) - 1/R' 


(4.13) 
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By use of the same deduction one can obtain the power lavs as 


f _ <J > 1/2 t " 1/2 and R _ ^ 


for low temperature case, and 


1/6 t 1/6 


f 


. ^ t' 3/5 and R . * 1/5 t 1/5 


for high temperature case. Similar but simplier analyses can be found 
elsewhere [ 21b). 

One may use Gibbs-Kelvin (Gibbs-Thomson) equation [46a], and write into 
the following form 

C(R) = C(<?°) exp[^/i(R)/kT] - C(<*> ) [ i+4/U(R)/kT] 

If higher order term(s) being added [46b], such as 
C(R) = C(Jo)[l + a_ 2 R _1 + a_ A R~^] 

with a ,=f( cT.v ), a =f(q,6, v ) [46b], In which CT js surface tension, q is 
— 1 m —4 m 


charge of the particle, approximated as 


2 

q^q^. =q , £ is dielectric constant 


of the medium, and v is molecular volume of the monomer as usual. However, 

m 

the surface tension is also a function of droplet size and is expected to 

decrease with decreasing in droplet size over wide range of circumstances 

[46cj, therefore the absolute value of a , and a , are governed by v , 

-1 — <* ~ m 

<T(R) , q, and 6 . Note that a , > 0, since 2<Tv > 0 [46d]; and a . < 0 

- 1 m ~- 

[46b] as long as q.c ; . < 0. For colloidal particles, always with residing 
charges, the above inverse power law in droplet radius is worthwhile to apply. 

IV. 6. System containing nonspherical sinks 


So far we only discussed systems with assumed spherical particles, 
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Odijk and Lekkerkerker recently reported [47] their work on theory of the 
isotropic-liquid crystal phase separation for a solution of bidisperse 
rodlike macromolecules. They used Gaussian trial functions that depend on 
adjustible parameters which were found by minizing the free energy of the 
system. Since the rods are never ideally hard particles, they are always 
flexible, therefore Odijk and Lekkerkerker took only the leading terms from 
the asymptotic expansion of the Onsager theory in studying the complicated 
influence of bidispersity. They found that the longer rods more highly 
ordered than the shorter ones [47]. 
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Table 2. Comparison of Selected Works on Solution of 
^C/^t * Dv 2 C + A, with A»0 and or A^O * 


Ref . 

B. C. and I. C. 

Solution, C(r,t), while A=0 or A/0 

Collins & 

C(R,t)«0 

C(r,t)“C o {l-(R/r)erfcl(r-R)(4Dt) _I/2 ]} 

Kimball 

c(<~, t)=c o 

Proposed B.C.: C(R)=( C/ r) , Y*=D/k 

K 

1949 [38] 

C(r,0)=C , r>R 

o ' 

R=sink radius 

D=dif fusion coeff., k=rate constant 
(See Appendix D.) 

Reiss & 

C(r ,0)=C 

(1) A=0 

La Mer 

C(x,t)=C o e~^ Bt 

C(r,t)=C e ^ Bt +C x/r£^ sink(r-x) 
o o n 

1950 

[42a] 

0C«r) =0 

r=h 

■/B=constant 
x=sink radius 
h=higher limit 
of boundary 

C(h,t)=<C( t)> 
(See Fig. 2 b) 

yB[exp(-PC 2 D t)-exp(-yBt)] *** 

(7B-* 2 D)tf 0 2 
n n 'n 

(2) A=-k(t)/D and C(x,t)=0 assumed 

C(r,t)=C x/r ^exp(-X 2 Dt)[sin< (r-x)]/p/ (3 
o n n n In 

+ higher order terms (x<r<h) 

i . _ - - 

Lif shitz 

c-c(t), c( 0 )=:„ 

j C -Cpo + ^/R, *(2 ff 7kT)v ffi C 

& Slvozov 

C =C at the 

| i= a5(r-C/<^r ) „=(&/R)(C - C ) 

1961 [9a] 

R ec , 

i * bounaarv 

C =Set. cone* 

R “critical 

c0 -j • 

radius at 
t=0 

p 

i 

‘ r=R R 

j =(£>/R)( - °^/R) , ^ “D/,5 , £ “density 

j =C-C^ “degree of supersaturation at t>0 

i 

1 /\ =C ~Cp==de£ree of supersaturat ion at t=0 
• ~-*c 0 

i 

dR/dt = £/F. ( /- - *Vr) 

! 
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Wagner | Oaverage cone. 

1961 i C '^surface cone. 

of particles 
[9b] I with radius=r 


C =saturated cone 
for particles 
with radius=r 


C =sat. sol'n of 
particles of r= o° 


C *C [1+2®V /prRT], >*>*] for undissociated 
r o a 

dr/dt=[kDv /(kr+D)](C -C) 
m r 

« -(Dv /r)(C '-C), for kr»D, dl-case 
m r 

= -kv (C '-C), for kr<<D, ik-case 
id r 

R=gas constant 
T=absolute temperature 


Felderhof 
& Deutch 
1976 
117] 


Marqusee 
& Ross ** *** 
1984 

[21b] | 


n( r)=0 , r=R 

=C, r= P° 

Here n(r) is 
equivalent 
to C(r). 


c (r,t) =C (R) 
K eq 

C( r , t ) =C 

r-*-t» o 


(1) A=0 , n(r) = n - n R/r 

o o 

J= f jds=4'*D Rn =k n , j=-D vn 
o o o o’ J o 

(2) A=-k Q ^<n>, ^<n>/^t=D Q V 2 <n>-k o ^<n> 

sink density 

3 

n(r)=n +e .r-n (R/r)-e ,r(R/r) 
o o o o 

(See Appendix B.) 

(1) A=0 , C(r,t)=C Q [ l-R/r]+C e ^(R)R/r 

J =410) R[C -C (R)] 

T o o eq J 

r*° 1 

(2) A=-4 ^/v^j R V[R,<C>]P(F.)dF 

V[R,C]=dR/dt=(D v /R)i<C>-C (F.)] 

O B SC 

C(r,t) see Appendix B. 

C(r,k) and C(q,A) see Appendix C. 


* A=sink/source term(s), B. C. =boundarv condition, 1. C. = ir.it ial condition. 

** Similar works in period of 1979-B3 have previously reported -;ic-2ibj. 

*** oC and B see definitions defined in [42a]-Eqs. 11 and 11. 
n | n 
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APPENDIX A 

EVALUATION OF THE GROWTH LAV 

Consider a spherical sink (the precipitate particle) of radius R, 

3 

(4fi/3)R /v E volume of the sink/molecular volume 
m 

= total number of molecules in one sink = n (per sink) 

Let I = dn/dt = rate of molecular adsorption or 
desorption on the sink surface 

= d/dt [ ( 4 1 T/3)R 3 /v ] 

m 

= (4-rr/v )R 2 dR/dt (A. 1) 

ZD 

3 

or equivalently, d[(4 ' r *’/3)R ]/dt = v dn/dt. 

m 

Note that, since, I = and J is the total flux of mass getting 
into a single isolated particle, i.e., 

J = - J ds.j = - 4~R"n.J (A. 2) 

s 

where ii is the out normal unit vector on the sink surface. The negative 
sign in Eq. (A. 2) taK.es care of the case for negative vC, in order to make 
j positive (see Appendix 5). By definition, 

j = - D ?C(r,t' = -D dC(r,t)/cT 
O o 

and since [ 2 i b . . 

C(r,t) = C [ I — R ■ r ] + C (R) R/r 
o ec 

is the solution of Lapiace-equation in the concentration field arround a 


43 



single sink, which satisfies boundary conditions 


C(r ’ t} r-R * C eq (R) and 
C(r,t) r=<)tf * C q , therefore 

j = -D [C - C (R))/R (A. 3) 

o o eq 

Substituting (A. 3) into (A. 2), one obtains, 

J = 4 7TRD [C -C (R)] (A. 4) 

T o o eq 

After equating (A.l) and (A. 4), one obtains for one single sink and consider 
t=0, 

dR/dt = (D v /R) [ C -C ( R) ] (A. 5) 

o m o eq 

For t > 0 and for a distribution of particles with number densioty, ^°, and 

probability size distribution function, P(R), we define an average 

concentration C (r,t) under the asumption of particle independence, 
av 

such that , 

C (r,t)/ t = D v 2 C (r,t) 
av o av 

- (£' r V7v ) | R 2 dR/d. P(R)dR (A. 6) 

' 11 'G 

The microscopic concentration fieic [ 2C , 21b] obeys the steady state 
diffusion equation (see Appendix B). 


Now, consider the rate of molecular absorption or desorption 
number-density of sinks n(?..r; or (R,t) as used by some authors, 
saturation { C— j , and the volume fraction ^ [14—16 ,18,21b] of 


K=J 

total 
the super 

the sinks , 


I = k(R) P = aUD R 

\ o o \ o 
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4^TD R(C -C (R) ] , as t«0, 
o o eq . ’ T 


(A. 7) 


I - k(R) 

- 4TTD R[C -C (R) ] , as t>0, (A.8) 

I ■= k(R) 

= 4 ^D q R [l+Rj4T5-^<R>nC av -C eq (R>] , as t>0,<j>>0 (A. 9) 

Here 4 TtD q R is responsible for the effect of the size of the single particle, 
[ 1+RJ4 T" ^ <R> ] is responsible for the effect of particle concentration, and 
[C -C ( R) ] is responsible for the effect of concentration field difference, 
the extent of supersaturation. 

The growth law now becomes: 

V(R,C )=dR/dt 
av 

= (D v /R) [ 1+R j 4> f <R>) | C -C ( R) ] ( A. 1 0) 

o m v \ av eq 

This is the growth lav and competition of sinks (precipitate particles) 
and under the influence of the concentration field gradient as well, now, 
the rate constant in (A. 9) becomes 

k(R) =4 "7TD R[ 1+R :TW<R>] ( A. 1 1 ; 

o ~ i 
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APPENDIX B 


CONTINUITY EQUATION AND ITS LAPLACE TRANSFORM 

As consider the sinks to be completely at rest, the monomer 
concentration, C(r,t), satisfies the diffusion equation 

^C(r,t)/^t = D Q V 2 C(r,t) (B. 1) 

in the part of the space occupied by the solution. The microscopic 
concentration field obeys the steady state diffusion equation 

v 2 C(r,t) = 0 (B.2) 

The so called Smoluchowski-boundary condition [24] requires C(r,t) to vanish 
at the surface of the sink, the problem has the obvious electrostatic 
analogue of finding the potential for a collection of ideally grounded metal 
spheres . 

According to Smoluchowski-theory [24] for the steady state reaction rate 
in a dilute system of spherical sinks, at sufficiently low sink density the 
competitiuon for solute molecules between sinks can be neglected and one 
needs consider only the effect of a single sink. One assumes that a monomer 
diffusing toward the sink is absorbed instantaneously upon contact and 
expresses this by requiring Cvr.t) to vanish at the surface of the sink. In 
a steady-state situation one tnerefore must solve Ec. (B.2). The boundary 
condition imposed on Ec . (B.2, is 

C(r,t) = C, at ‘r-r. d= ?. (B.3) 

i 

where R is the sink radius (see Fig. 2). The initial condition is 
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C(r,t) - C Q (r) at t - 0 (B.4) 

For uniform concentration C q at infinity, the solution is 

C(r) = C q - (C o R/r) (B.5) 

For convenience, one may choose the coordinate origin located at the center 
of the sink, so that from r=0 to r*infinity, there are numerous vectors with 
r=R located on the surface of the sink. In this manner the meaning of k Q 
is self-evident, since 

k Q » 4 Tr D Q R implying k Q proportional to R (see below) 

The molecular flux encountering the sink surface is 


j = -D^vC , by operating v on Eq. (B.5) one obtains 

vC * C R/r 2 
o 


J = j j.ds = the total flux absorbed per second per sink 
s 

after integrating over the surface of the sink 


i j(D C R/r“ r~ sin© d© d<f-- 

J. ■/' o o ' 1 

J c c 


= 47TD RC = k C 

CO CO 


Or.- a macroscopic scale rhe equation for the average monomer concentration 

is C = <C> , therefore the continuitv eouation becomes 
av 

<0/>t = D r~<C> - k f <C> (E.c.- 

o o \ 


with the rate constant k 

o 


4"wD R and sink number density 
o 


r- 


Since C(r,t) in a concentrated solution varies rapidly, therefore it is 
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no longer true that the concentration tends to a constant at a large distance 
from a single sink. The magnitude of the effect of gradients may be 
estimated by solving Eq. (B.2) with the condition 

C(r) = C q + e Q * r as r approaches infinity [17] 

The solution is: 

.. C(r) - C + e .r - C (R/r) -e .r(R 3 /r 3 ) (B.7) 

o o o o 

The particle flux into the sink is still J = AttD^R, so that it appears that 
concentration gradients have no effect. However, the dipolar disturbance 
(the last term in Eq. (B.7)) makes itself felt at large distances and, 
when summed over a large number of sinks, contributions significantly to the 
average local particle density. 

If replace Eq. (B.l) by a microscopic equation valid for all space with 
added sink terms representing the monopole and dipole contributions [17], 

•> K 

<?C(r k t)/it = D [ v t C( r , t ) - q.$(r-r ) 

j = l 3 3 

N 

+ p .v$(r-r )j (B.S: 

j = l J J 

where C(r,t)=0 a: r-r . <=F. (inside the sink or at the sink surface.', Bv 

3 

use oi (B.S) one may take average as follows: 

LKS: <^C(r,t).Vt> = <fC ,'J t = c t 

£ v- 

? N . N 

RHS : <D "C(r,t)>=D v“<C>-4»D <Z- c . £(r-r . )>+4 D < p . . vr ( r-r . )> 

O O O '1 6 j 0 J ° J 
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After equating the two sides, one obtains: 
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<^<C>/c)t - D v 2 <C> - 4ITD Q(r,t)+47ID P(r,t) 


v.(D(^)v<C>] - k(p f <C> 


(B.9) 


with k(j')f<C>“4in) o Q and D(pv<C>=D ( V<C>+4'nD o P [17]. Eq. (B.9) accounts 
for the competition for monomers between sinks, as implied by k(^). It 
also shows, by D(^), the diffusion motion will be affected by sink 
concentration ^ . 

If choose only the monopole sink term which only approximates the 

boundary condition but not affects the first order correction. The steady 

2 

state solution of D q v C(r,t)= (r)+ q^ (r-r^) has the formal solution [21b], 

C (r,t)= G o*4U> +lG o (r-r.)q. 

Here the symbol G * denotes the convolution, as an integral operator, 


G o *f( r )=J dr ' G 0 ( r - r ') ( r ') 


and G (r)=-l/41TD r is the free diffusion propagator. The Green's function 
o o 

2 

for the operator D q v is, according to Creen's Theorem [23a, p.242]. 


;( r-r ' ) = 


- 1 / ( 4TB Ir-r'i ) 

o . 


Actually, "Smoluchowski" boundary condition, CfK.t) = 0, at the surface 
of the particle is not self-evidently valid [25a]., so that the "radiation" 
boundary condition [25b; may be employed. Wiiecsk: ar,c Fixman [25c] compared 
these two boundary conditions and made them equa^ as set 6- =Q in C(R+f.. 


The continuity equation (B.l) with boundary condition and initial 
condition (B.3) and (B.4) respectively is more convenient to work with the 
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Laplace representation obtained by Laplace transformation of Eq. (B. 1) as 
follows : 

L f(t)=Laplace transform of f ( t )= J f ( t )e"^‘ t dt=f^(-^) 

2 

L {^C(r,t)/<*t * D q v C(r t t)}=Laplace transform of Eq. (C.l) 

. ao .OO 

LHS : L <)C(r ,t)/J t =f^C(r,t)/^t e^dt = C(r,t)e" Xt + ( Xc( r , t )~ Xt dt 

O 'c 'o 

= -C o (r) +Xc(r,X.) 

OO 

RHS : L D v 2 C(r,t) = D v 2 [ C(r,t) e _Xt dt = D v 2 C(r,X.) 
o o o — 

The Laplace representation for (B.l), with boundary condition (B.3) and 
initial condition (B.A), now becomes: 


2 

jC(r,*-) - C Q (r) = D o v £(r,X.), \ is the Laplace transform variable, 


C( 


r-r . 
1 


<=R, M = 0, i=l ,2, K. 


This difussion problem can be viewed as a case of multiple scattering, 
and the above equations are equivalent to a set of simultaneous equations by 
statistical averaging and mean field approximation [17,18] or by Fourier- 
Lapiace transform tc obtain a configurational average Green's function [19]. 
The latter method has been usee by many aurthors [ 2 C , 2 1 ] . A brief discus- 
sion is shown ir. Appendix C. 


52 



ORIGINAL PAGE IS 
0E POOR QUALITY 


APPENDIX C 

MULTI SCATTERING REPRESENTATION AND ITS FOURIER TRANSFORM 


As shown in Appendix B, the Laplace transformed representation for the 

2 

partial differential equation, ^ C(r,t)/^t = D q V C(r,t), is an ordinary 
differential equation in the form 

C(r,X) - C (r) = D v 2 C(r,X) (C.l) 

— o o — 

where X is the Laplace transformation variable amd £=C Q (r) at t=0 and 
C( | r-rj <=R^ , ^')=0 are initial and boundary conditions respectively. 

Where i=l,2 N for the radii of the N traps (sinks) located at {r^} 

respectively. Many authors use "R” or "a" for all Jr-r^j > just for 
convenience. Rearranging (C.l) we obtain: 

(X- D o v 2 )C(r,A) = C Q (r) 

C(r,\) = (X- d v 2 ) *C (r) 

— o o 

= G (r)C (r) in absence of sinks, 
o o 

If more than one sinks are considered, one should ado the sink terms. In 
this case the solution becomes: 


X' = 


r ! 


N 

i = 1 


( r ) t . C , ( r ,X) 


where t, is me t-operator for one sink located at and _C(r,X) is the 

Laplace transform of the concentration in the system which the sink at 

is eliminated. G (r) is the Green's function in the absence of traps, 
o 

According to Green's theorem in handling boundary conditions in 
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electrostatic or analogy [26], a simple application of the divergence theorem 
leads to the solution of D Q V 2 C(r) = f(r) being C( r )= J G( r ,r' )f ( r' )dr' , 
provided G(r,r')=[D Q v ] £(r-r') is imposed, where £(r-r') is the 

Dirac delta function. The solution of 


N 


V - r,X ) = 4(r)+Z q i ^(r-r i ) is 


f N 

C(r, M= ] dr'G(r ,r')[ 4 (r')+ $(r'-r )] 

y i=l 


The physical meaning is that (Xr,^-) can be obtained by integration with 

2 — 1 

an inverse operator of [A -D q v ] . The inverse operator is equivalent 

to Green's function in such a way 


C(r,K) = i G(r,r')(4(r')+ X q £ ( r' -r )] dr 


i=l 


JlA-c^r 2 ] 1 [ < f(r')+ X q ± ~ r i ^ ^ d r ' (C.2) 

i= 1 


In (C.2), only monopole terms are considered. Replacing (A-D_v j * d\ 
G Q (r) and introducing another t-operator, one has 


Ci-.K - G (r)C (r) t Y, G (r)t.C.(r,A) 
— o o : — . o i—i 


..Is, 


= G ( r ) C (r) -f Z_. G (r)t.CL(r, ) 
C* C 1 , . C "I ; 




where C__.(r,'*~: is the Laplace transform of the concentration in the system 

from which the sink at r. is eliminated, and t. is the t operator for one 

I I 

sink located at r^. 
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By substituting (C.3b) into (C. 3a) one obtains the multiscattering 
[ 18-20 ,21b, d] representations , where is an integral operator [23e], 

N 

C(r,X) = G Q (r)C o (r) + ZL G o (r)t i {G o ( r)C Q (r) 


N 

+ 2Z G o (r)t.C.(r,X)) (C. 3c) 

J ^ 9 ^ ^ 


SAc (r '’ X) = jv r ’ r '* X;r k^ (r ' X)dr ' 

t.C.( r "A) = T j (r',r^;r j )C j (r" A)dr" 


To solve (C.3a) it is convenient to use Fourier transform technique. The 
Fourier transform of a function, say f(r), is defined as 

DC 

F f (r ) =J f(r)e iq * T dr = f(q) 

— Oc 

therefore the Fourier transform of Eq. (C.3a) is 


LRS: F C(r,M - ( C(r, )e 1C; ‘ r dr = C(q, x ) 

J -po 

A/ 

RHS : F [G (r)C (r) + Z. G t. C .(r, x )] = (I) + (II), wher. 

G O 0 1 — 1 


(I) = \ G ( r) C (r)e iqr dr 
o o 


i r 2, icr, , 

= i.i * r jcr/ : ( K -0 v je j : 
o o 


• | /« ' . f i r f • \ iG* • \ ^ ^ 

= ‘ {C irjcr. [\-ii (lq) e - j;, \iq> = -c 

• G C 


■ S n r 

SC, (1) = ;C (r)iA-rZ' C~; *fe '‘dr 


page 13 


" _ .IIS! 

QP POOR Q t v 


,ur 


= • ■ 


, . -icr 

t „ •. r ; e ■ g r 


= G ( q)C (q), where G (q) = iA+D c‘] *, similarly 
OG O G 
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= L G (q)T (q)c (q» 

i= l 

After equating both sides, one obtains: 


C(q, x ) = G (q)C (q) + £ G (q)T (q)C (q» x ) 

i=l 

and T i (q) = (3 ) 3 j dq'T^ q ,q' , X jr^Ptq ,q' ) 


(C.4) 


T^(q) is an integral operator, in which P(q,q') is a permutation operator 
that changes the index q of the operand to q' , and T.(q,q', ; r ^ ) is a 

convolution integral operator, in Fourier transform, 


^(q >q' , x ;r 1 ) = jdrjdr'e l(qr q r ^T.(r ,r' , * ;r i 


) 


. ^ ^ 

Substituting C(q , x )=G o ( q )C q ( q ) + G q ( q )1\ (q)C^ ( q ,^) 
into (C.4), one obtains 


~ N 

C( q ,> ) =G ( q ) C ( q ) + I_ G ( q ) T ( q ) { G ( q ) C ( q ) 

O C G i O O 

1 = ‘ 


G ( c, )T . ( c )C . ( c , x ) } 

O' j 's' 


(C.:- 


ir. suer, a way one may substitute 


C_. (q. ; S'=G, (c JC .(c)* G (■:)!. (q)C. (c into (C. 5) , and sc on. 

“i C' G O k K. 


K = i r =r'; 
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APPENDIX D 

SOLUTION OF FICK'S 2nd LAW, EXAMPLE 


The solution of 2nd order partial differential equation, Fick's 
second law of diffusion, 


^ C/2t *= D v 2 C 


(D.l) 


chosen by Collins and Kimball [38a] with the boundary condition CCR.O^O on 

the surface of a 6ink of radius R (see Table 2), is similar to the solution 

1 /2 

of linear flow of heat [41]. For the normalization constant (4 7TD) , 
the particular solution should be like this: 


n 2 

C(x,t) = t exp(-x / 4Dt) 


2 2 1/2 

since by change of variable, x /(4Dt)=u , du=dx/(4Dt) , 


(D. 2a) 


-pc 


( C(x,c)dx = 2 j C(x,t)dx (since C(x,t) in (D.2a) is an even function), 
* J c 

= 2C7TD)*' /2 f exp(-u 2 ) du 
•x 

= 2( TT D) 1/2 = (4TD) 1/2 ( D . 2 b ) 


1/2 

(4 T D) is the normalization constant or the chosen function, the error 
function solution. Now, consider one spherical sink of racius R and locate: 
at the origin, then the error-function-solution takes the following form 


C(r,t) = C {] 
o 


- ( R ' r / 


;rfe [ ( r-F.) .'(4Dt ' 1 j i 


( D . 3 ) 


The solution satisfies the bounaarv conditior 


C( r , t) = 0, for 0 <= r <= R 
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so that the concentration field of the monomer is only defined in the range 

of r-R to . In other words, the concentration of monomer is not defined 

in the region from r=0 to R, i.e. , inside the precipitate-particle of 
3 

volume A'hR / 3 , since inside this volume monomers not exist. 

The total flux J, which is J^j.ds over the entire surface of the 
assumed spherical precipitate particle, was defined by collins and Kimball 
as £ , i.e., 

? = J = j.ds (j=-DvC, ds is negative to j in this case) 

= 4H R 2 D^CAM 

r~ k 

= 4~R 2 D^/^r C {l-(R/r) erf c[ ( r-R) /( 4Dt ) ^ 2 ] } 

o r = R 

= 4-rrR 2 DC {(1/R) erfc[(r-R)/(4Dt) 1/2 ]} D 
o r-K 

+ ((R/r) erf [ ( r-R)/( 4Dt ) 1 /2 J } __ (D.4) 

r — k 


To obtain Eq. (D.4), the following Leibnitz-rule [44] should have been used, 


d/dt 


b(t) 


/ a(t) 


f(x,t)a>: = f [b( t) ,t]b'(t)-f [a(t) ,t]a'(t) 


f ( x , t ) ic - 1 j dx 


a\t) 


Besides, by definition. 


erf (>:) = ['- . ~ ?T ; " j exp )c>:' 

erf c (x) = I - eric x 


one can see that 



erf (0) * 0, erfc (0) * 1 - erf (0) * 1 
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erf ( p ° ) = 1, erfc (• > °) = 1 - erf (p° ) = 0 

1 / 2 

Since r-R=0 if r=R, so that erf (r-R)/(4Dt) = 0 as r=R, therefore 

we obtain a very neat fora for the total flux as follows: 

= 4TTDRC o ll + R/(7TDt) 1/2 ] (D.5) 

Moreover, Collins and Kimball defined sk (=J) as the time derivative of the 
total amount of the diffusion species, from the amount presented at t=0. i.e., 

(=J) = dM/dt 

We can see from Eq. (D.3), the total amount of the precipitate is just equal 

to the integration of the second term on the right hand side of Eq. (D.3) 

over the entire region of the precipitate-particle sizes, for spherical-size 

2 

assumption, the volume element is dv'=4Tf'r dr, i.e., 
f ^ 

M = C ) 4*r~{(R/r) erf c [ { r-R) / ( 4Dt ) 1 ] } dr 
° 'R 

= C 4 t r r p, erfc[(r-R)/(4Dt)^ 2 ]dr (D.6) 

'R 


which is the total amount of precipitate in all sizes of particles at time 
t>C>. By taking time derivative of Eq. (D.6), and by use c: Leibniz Rule 
again, one ootains 


dM dr = 


d/ct-C 4T r R e 


rfc!(r-R)/(4Dt) 1/ -jdr 


f 




C 4-^r; r(d/dt){2/( ~ r T) il M 

° 'R \ r-R) / ( 4Dt ) 


extv )ax }dr 
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C 4'tnU4Di<)~ 1/2 t‘ 3/2 
o 


( 


r(r-R)exp[-(r-R; /(4Dt) ]dr 


2 - 1/2 

then, by change of variable: (r-R) *=u, dr=(4u) du, and change of 
integrating limits from r=(R, °°) to u=(0, e ’°), then 


,tx> 

£ = dM/dt = 27TRC (4DT0* 1/2 t~ 3/2 { I u 1/2 exp[-u/(4Dt) ]du 

0 Jo 


+ R j exp[ -u/( 4Dt ) ]du) 

J 0 


(D.7) 


Eq. (D.7) can be evaluated by use of the following integrating formula, 

OC- 

( i /2 _ 1/2 

J x exp(-ax)dx = ( n /a) l/(2a) 

c 

provided we set the constant "a" equal to (4Dt) *. The integrated form 
of Eq. (D.7) is. 


±= dM/dt = 2-KRC (4Dff) i/2 t 3/2 [ 2Dt ( 4Dt TT ) 1 /2 +4RDt ] 
o 

= 4T7RC DlD~t)~ 1/2 l(D TTt) 1 ' 2 + R] 

G 

= 477JJC [1 + R/(Drt) 1/,: ] (D.8) 

G 

Eq . (D.8) is just me total flux. £. , represented by Ec . ID. 5). From th 
above argument one can see tnat the significance of the solution of the 
second order cirfcrsniiai coucLi on of Picks lav of diffusion chosen by 
Collins and Kimball :38a':. 
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a) a_,R 1 can be derived, from the principle called equilibrium 

configuration according to Gibbs-Wulff (see (9bj-p.582). By this 

principle the exponential function in Gi bbs-Thomson equation car. 
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vr.ere (f is the inter- facial tens ion , e> Ole V=t)( 47TR“ )/<)( 47TP /3)=2/F. 
is the quotient of surface area changes and volume changes during 


ne growth or the crvstal, k is the Boltzmann constant and r. is 
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racius of the assumed spherical crystal grain. 
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FIGURE CAPTIONS 


Fig. 1 
Fig. 2 

Fig. 3 

Fig. 4 
Fig. 5 
Fig. 6 

Fig. 1 
Fig. 8 

Fig. S 


. Cutoff - pair potential, which vanishes as r >= b. 

. (a) Two dimentional schematic representation of circular particles 
in Dirichlet region. 

(b) Reiss & LaMer's model [42a], h-x = thickness of the diffusion 
region (see Table 2). 


Schematic representation of N(t) = n(R,t)dR. See Section III.l, 
and Kahlweit and Teubner's 1976 and 1980 works [22aii, 22b]. 


a )( z , ) vs z, ik-case (see Table 1). 


4b’( z , (j- ) vs z, dl-case (see Table 1). 
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Cutoff size distribution, z vs 
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Reduced coarsening rate vs equilibrium volume fraction^ vs 

*■» n 

^ ^ po ) = d a (, ) ■ i d a ( 'j ) ) { 2 1 b j • 


Amplitude for the density of particles f 0 'f‘ 


pc ) vs 


ftc , 


3 ver< 


radius a (4^ ) vs ^ , and the supersaturation decay function 


vf 


u v t> ;* v - 


ln which a ( 0 ) = ( 4 / 9 ) 1 / ^ , Cr(0) = (9.-' :/ ', 


I** 


and = 1.99 are usee [21bj. The t^ /J , t i '~, and 


t J factors are taken as II + 0.81 5(-f<® ) 
1/2 - 


i / : 


, [ i - 0 . 8 1 3 ( ^ •* ) ‘ “ j , 


and { 1-2. 44 5 (^ f 


i resoect ivel v. 


69 






TAt. *t.cL*^tr 


(£i) £ - *//>» #*J tcmt^L. Xc/-S- Ma Ci 

?t.jb'l}£rt 7i-/tcr>) Of- C<' TCu-f C*T_ 

$ /« 2?’' rict /*-£- 




(_6) /?e/jx ^ /y £r i 

TvodeJ <*-] , 

h - X- — '’CS’S 

<y- tu </•//« 

-<-*-« /<rrs_ C*jut- Z06/C *Z )_ 





u)( 25) ys 2? i k — case 




••*7 = v'Z 

V 'C 

°7 =/ 



$<*!*■ cC y ^ 



c*J 








ORIGINAL PAGE *13 
OF R^OR QUALITY 


Amplitude for the density of parti 
ties us phi, auerage radius us phi 

Oc Sigma O ‘'firsat. -fu>/cX'o7^ } pH j . 






